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Abstract. The SUE machine is used to extend by a factor of 1000 the time-scale of 
previous studies of the aging, out-of-cquilibrium dynamics of the Edwards-Anderson 
model with binary couplings, on large lattices (L = 60). The correlation function, 
C(t+t w ,t w ), t w being the time elapsed under a quench from high-temperature, follows 
nicely a slightly-modified power law for t > t w . Very tiny (logarithmic), yet clearly 
detectable deviations from the full-aging t/t w scaling can be observed. Furthermore, 
the t < t w data shows clear indications of the presence of more than one time-sector 
in the aging dynamics. Similar results are found in four-dimensions, but a rather 
different behaviour is obtained in the infinite-dimensional z = 6 Viana-Bray model. 
Most surprisingly, our results in infinite dimensions seem incompatible with dynamical 
ultramctricity. A detailed study of the link correlation function is presented, suggesting 
that its aging-properties are the same as for the spin correlation-function. 



Submitted to: J. Phys. A: Math. Gen. 

PACS numbers: 75.10.Nr,75.40.Gb,75.40.Mg 



Aging in Spin Glasses in three, four and infinite dimensions 



2 



1. Introduction 

Spin-glasses [U El EH E| were discovered to age even on human time-scales some 
twenty years ago [3]. Aging is nicely demonstrated, for instance, in measures of the 
thermoremanent magnetization (see e.g. jH]): in the presence of a magnetic field, cool 
an spin-glass from room temperature to the working temperature, T, below its glass- 
temperature; hold the magnetic field for a while (the time elapsed will be called t w 
hereafter), then switch-off the field and record the time-decay of the magnetization 
M tw (t). Not only this decay is very slow but, even for the longest t w tried up to 
now, M tvj (t) strongly depends on t w (the larger t w is, the slower decays M tw (t)). It 
has slowly become clear that the important information coming-out from experiments 
in spin-glasses regards dynamic out-equilibrium effects, such as this one or the more 
sophisticated memory and rejuvenation effects 0E]. Although there has been a burst 
of theoretical activity in out-of-equilibrium dynamics [HI 13 UH] , it is still not clear:): how 
out-equilibrium experiments will help us to choose among the conflicting theoretical 
views on the nature of spin-glasses: Replica-Symmetry Breaking :3j (RSB), the droplets 
picture ^3], and the intermediate TNT picture ^H]. The situation is farther complicated 
by the fact that detailed theoretical predictions (to be confronted with experiments) can 
be extracted from models only through Monte Carlo simulationsjn]. It is worth recalling 
that numerical results on out-of-equilibrium dynamics have cast some doubts even on 
the usefulness of the Edwards- Anderson model to describe physical spin-glasses [TH] (see 
however for some reassuring results). 

It is thus clear that one needs to address quantitatively the time-decay of M tm (t) 
or equivalently, given the Fluctuation-Dissipation Theorem jH], the time-dependent 
correlation function§ in the absence of a magnetic field: 

C(t,t w ) = j-Y,(o- i (t + t w )a l (t)) (1) 

i 

Now, it seems to be a fact of general validity in out-equilibrium dynamics |S that 
C(t,t w ) behaves differently in different time-sectors. Loosing generality|| for the sake of 
clarity, this amounts to say that it can be decomposed as 

C(t, t w ) = £ /, ( (t + tw) 1 1 ^~^ i ) . (2) 

Here, fi are smooth, decreasing functions that tend to zero at infinity, and such that /;(0) 
is of order one. It follows that if t <§C t% the i — th time sector contributes the constant 
value /i(0) to C(t, t w ), while for t 3> t^ it contributes nothing. In other words, the i— th 
time sector is active only for t ~ tf£ (notice that the different time sectors get neatly 
separated only in the limit of very large t w ). Not much is known about the exponents 

| A very encouraging experiment measuring the violation factor of the fluctuation-dissipation 
theorem ^2 ^] has been recently reported. 

§ Actually, in the aging regime they are related by a very smooth function ^2 Q2| • 

|| The here presented formulation cannot describe logarithmic domain-growth, for instance. See in |HJ 

the general framework. 
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fii defining the different time sectors. With this popular parametrization (HUH], one has 
< fa < 1. For the simple case of the coarsening-dynamics (domain-growth) of an 
ordered ferromagnet |18j . only two time sectors are needed for a complete description: 
Hi = describing the stationary, t w independent dynamics found at small t, and /i 2 = 1 
describing the full- aging situation where the correlation-function depends on the ratio 
t/t w . Also the spin-glass dynamics has been experimentally claimed^ to be ruled only 
by two time sectors: fa — and fa = 0.97. The second time sector is slightly but 
clearly different from the full-aging t/t w behavior (fi = 1) and is thus named sub-aging. 
However, a very recent experiment |1 9 j seems to indicate that the sub-aging behavior 
is just an artifact of the finite-time needed to cool the system down to the working 
temperature (a limitation not suffered of in numerical simulations). Using their fastest 
cooling protocol, Rodriguez et al. [TH] have found a clear full-aging behavior /i 2 = 0.999. 
Furthermore, the role of the stationary time-sector (fa = 0) to describe the data is far 
less critical than previously |)j thought. It is also worth mentioning the recent numerical 
results in 3 and 4 dimensions by Berthier and Bouchaud |T2|, who found superaging for 
infinite cooling rates, turning to subaging for finite cooling rates. At this point, we wish 
to make two comments: 

• The presence of more than two time-sectors seems to be a crucial requirement 
for the validity of the dynamic version of the usual ultrametric Replica-Symmetry 
Breaking description of spin-glasses with an infinite number of replica symmetry 
breakings. 

• If the largest [i exponent is 1, as implied by this popular parametrization |UJ |S] , for 
the very large t w achieved in experiments (1 second means roughly t w = 10 12 ), it is 
quite possible that all the faster sectors have already died-out since t% <C t w ~ t. A 
full-aging ansatz could pretty well describe the data, specially if the second-largest 
/i exponent is significantly smaller than one. In this respect, a numerical simulation 
could have a better chance of observing the different time sectors. 

To conclude this introduction, let us recall the main results of previous intensive 
numerical studies of aging dynamics jH] (for a more recent extensive numerical study of 
the aging dynamics with a different focus see [IE])- Indeed, it was found [5] that the 
correlation function at short times (t <C t w ) could be described as a power law with a 
temperature dependent exponent (C(t,t w ) ~ t~ x ^). The exponent x(T) was found to 
be fairly small (x(T = 0.5T C ) = 0.015) On the other hand, at long times a power-law 
decay with a different exponent, A(T, t w ), was observed (C(t,t w ) = t~ A ( T '*») for t 3> t w ). 
Yet a cross-over functional form was proposed [20] (see also Ref. |3]), implying that 
\(T,t w ) saturates to A(T) for large t w : 



The cross-over function was proposed to be decreasing, smooth, and to have the 
asymptotic behaviors $(y) ~ constant for small y and $(?/) ~ y-K T )+ x i T ) fo r large 
y. However, in the spanned time scales P (t w < 10 5 , t < 10 6 ), A(T) was actually found 




(3) 
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to significantly depend on t w . Yet, if one assumes that \(T,t w ) saturates to A(T) for 
large t w , so that Eq.(j3J) could hold, one would speak of a dynamics with effectively three 
time sectors: ll\ = 0, Li2 = 1 — x(T)/X(T), and LI3 = 1 (to derive this, assume that 
$(?/) is an analytical function, such as (1 + y)~ A ( T ) +:c ( T )). However, since x(T)/A(T) was 
found jH] to be 0.1 or smaller, it would not be easy to separate the Li2 and ^3 sectors. 
On the other hand, it is clear that Eq.Q describes a non-ultrametric dynamics. 

In this paper, by extensive Monte Carlo simulations using among others the 
dedicated SUE machine [2HE21I2H!, we wm show that there are small yet quite significant 
deviations to Eq.Q in finite dimensions. Surprisingly enough, our results in infinite 
dimensions are rather different, and suggest that dynamical ultrametricity could not 
hold in (some) statically ultrametric systems. As the reader will notice, heavy use of 
data fitting will be made in the following. Unfortunately there is no precise theory that 
tell us which should be the precise functional form of C(t + t w ,t w ), so it is difficult 
to justify theoretically many of the fits. Also the very good power like behaviour we 
found cannot be recovered from analytic computations. Here we are doing some kind of 
exploratory work, trying to guess which could be a reasonable form that well represents 
the data. This may be useful for a number of reasons: If we can do good fits at given 
vale of t w it is rather more convenient to look to the dependence of the fits parameter 
on t w than to the data themselves (sometimes scaling plots may be misleading), good 
fits may evidentiate some behaviour that could be eventually analytically derived, fits 
can also used to extrapolate data to large times. Finally, notice that as an outcome 
of this strategy, some of the fits that we are doing could also be useful in analyzing 
experimental data. 

The layout of the rest of this paper is as follows. In section |2] we describe 
our simulations. In section |3] we concentrate on the spin-spin correlation function, 
presenting a new parametrization of the function C(t, t w ), and discussing the possibility 
of numerically studying the existence of more than two time-sectors. In section H] we 
focus on the aging behavior of the link-overlap and the link-correlation function (defined 
in section EJ). We shall conclude that even in the limit of infinite waiting time, the link- 
correlation function ages. Finally, we present our conclusions in section 

2. The simulation 

We have studied the three dimensional Ising spin glass defined on a cubic lattice 
(L x L x L) with helicoidal boundary conditions^. The Hamiltonian is 



The large time scale and lattice sizes simulated have been possible due to the use of 
a dedicated computer, SUE (Spin Update Engine, Universidad de Zaragoza) based on 

% Let (x,y,z) be the lattice coordinates of spin number i = x + yL + zL 2 , then the coordinates 
of the three nearest neighbors (in the three positive directions) are given by: i\ = (i + l)modL 3 , 
i 2 = (i + L) mod L 3 and i 3 = (i + L 2 ) mod L 3 . 



n 




(4) 



<i,j> 
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programmable components, and achieving an update speed of of 0.22 nanoseconds per 
spin. Details about the machine can be found in Refs. |2"T) 12*2*1 ^1- 

The volume of the system is V = L 3 , Oi are Ising variables, (uncorrelated 
quenched disorder) are ±1 with equal probability, and the sum is extended to all pairs 
of nearest neighbors. The choice of helicoidal boundary conditions is mandatory (for 
us) because the hardware of the SUE machine has been optimized for them. During the 
simulation we have measured the following quantities: 

1 E (o-i(t + t w )o-i(t)) , 



C(t, t w ) 

Clink (t, 

q{t w ) 



N 



D 



DN 



E E + O + tw) 0-i(t w )0- i+fl (t w )) 



i 



D 



(2), 



.(2) 



DN 



(5) 
(6) 
(7) 
(8) 



in the above equations, D is the spatial dimension, stands for the nearest neighbor in 
the /i direction, while the superscript (1) and (2) refer to the replica index (real replicas: 
pair of systems evolving independently with the same couplings Jij). The notation 
(. . .) refers both to average over thermal histories (random- numbers) and over disorder 
realizations. Let us just recall that since the starting configurations are random, we 
have explored the so-called q = sector [""], in which one expects dynamical correlation 
functions to be self-averaging. 

Given the unique features of the SUE machine we have preferred to use it for 
very long runs in a rather small number of samples. The lattice sizes studied have 
been L = 20, L = 30 and L = 60. We have considered three values of the inverse 
temperature: (3 = 1.25,1.67 and 2.0 (hereafter to be referred to as T = 0.8,0.6 and 
T = 0.5 respectively). The critical temperature for this model is T c = 1.14(1), thus 
the selected temperatures are 0.7T C ,0.53T C and 0.44T C respectively. Given the very 
slow growth of the spin-glass coherence length (see e.g. jJ3"211) one should not expect 
noticeable finite-size effects even for the L = 20 lattice. However, as it is well known, 
sample-to-sample fluctuations in C(t,t w ) decrease fastly with growing system sizes. 

We can sum up the details of the SUE simulations in table [2 The SUE time- 
step corresponds to 8192 full-lattice sequential heat-bath updates. For C(t,t w ) and 
Ciink(Mio) we have selected the t and t w values in a logarithmic scale. These values 



corresponds actually to 8192 x 2 



0.25i 



The number of simulated samples has 



been 16 for the L = 60 systems and 32 for the L = 20, 30 ones. Here, we will only present 
the results for L = 60, since the data for L = 20 and L = 30 are fully compatible with 
them, but far noisier. Yet, much more accurate data have been obtained for moderate 
t and t w from simulations on a PC (see below) of an L = 33 system. 

To have some data at times shorter than SUE step 8192, we have simulated (using 
heat-bath) on a personal computer 80 samples with L = 60 (t w < 6877 , t < 65519). 
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L 


T 


Number of iterations 


Number of samples 


60 


0.5 


10 9 


16 


60 


0.6 


8 x 10 8 


16 


60 


0.8 


5 x 10 8 


16 


30 


0.5 


5 x 10 s 


32 


30 


0.8 


5.5 x 10 s 


32 


20 


0.5 


2 x 10 9 


32 



Table 1. SUE simulations details 



t =8192 




10° 10 1 10 2 10 3 10 4 10 5 



1+t/t 



Figure 1. Double logarithmic plot of the correlation- function obtained with SUE 
using a heat-bath algorithm, versus 1 + t/t w , for t w — 8192 and T = 0.5,0.6 and 0.8. 
A single power law fits satisfactorily the data for t > t w at all-three temperatures. 



In addition, we have performed Metropolis simulations on personal computers for 
the same model in 3D [L = 33, T = 0.8 = 0.7T C , t w < 2 18 , t < 2 22 ), 4D [L = 37, 
T = 1.2 = 0.6T C , t w < 2 18 , t < 2 22 ) and in the infinite-dimensional Viana-Bray model 
(z = 6, N = 5 x 10 6 , T = 0.8 = 0.34T C , t w < 2 13 , t < 2 17 and N = 10 6 T = 0.8, 
t w < 2 19 , t < 2 16 ). Even if the times were shorter than in SUE, the number of simulated 
samples has been much larger. The typical statistical error in correlation functions (as 
calculated from the sample-to-sample fluctuations) obtained in the PC was 20 times 
smaller than the statistical errors of the SUE results. 
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Figure 3. As in FigEl for t w = 39813120. 



3. Aging dynamics 

We have found in finite dimensions (D = 3 and D — 4) that the correlation function 
C(t, t w ) can be nicely fitted for t > t w as 

/ f \ -l/a(tte) 

C(t, t w ) = A(t w ) ( 1 + — J , t > t w . (9) 




Figure 4. Top: exponent a(t w ) defined in Eq.lQ versus t w for the heat-bath dynamics 
of the Edwards- Anderson model in D = 3, at 0.7T C . Data for t w < 8192 were obtained 
on a PC, while for larger t w correlation-functions from SUE were used. Bottom: same 
as in top panel for the prefactor A(t w ), defined in Eq.@. The solid line was obtained 
fitting for t w > 8192. 



As it can be seen in Figs. II 121 and El for a wide range of t/t w and t w and all three 
temperature a simple power-law decay seems enough to describe the data, although the 
coefficients A{t w ) and a(t w ) clearly depends on temperature (notice that a{t w ) is just 
the inverse of Rieger et al. [H] X(T,t w ) exponent). The behavior in infinite dimensions 
(z = 6 Viana-Bray model) is rather different and will be discussed at the end of this 
section. 

It is clear that the waiting-time dependence of the prefactor A(t w ) and the exponent 
a(t w ) is of utmost importance. Should A{t w ) and a{t w ) tend to constant non vanishing 
values, dynamic ultrametricity [S] would not hold, implying that the usual dynamic 
formulation of the ultrametric approach to continuous Replica Symmetry Breaking 
should be modified. Furthermore, we have been unable of finding a divergence law 
for a(t w ) compatible with dynamic ultrametricity, if A(t w ) has a non- vanishing limiting 
value for large t w . Thus, we tend to believe that dynamical ultrametricity implies that 
A(t w ) should vanish in the large t w limit. 

To obtain the coefficients A(t w ) and a(t w ) we have fitted the correlation functions 
C(t,t w ) to the functional form in Eq.(JHj). Yet, although Figs. IH2l andlo1 suggest a pure 
power-law behavior, we have found some dependence of a{t w ) on the fitting-window, 
particularly for T = 0.5. Specifically, at T = 0.5 a(t w = 8192) grows a 10% if the fit 
is performed for t w < t < l^t w as compared to the fit in the window t w < t < 10t to . 
Not taking care of this could be dangerous since, obviously, the larger t w the shorter 
the achievable t/t w . Thus, in order to isolate the ^-dependence we have restricted 



Aging in Spin Glasses in three, four and infinite dimensions 9 

T = 0.6 




10 2 10 4 10 6 10 8 



Figure 5. As in FigHfor 0.53T C . 

ourselves to the fitting range t w < t < 10t w . On the other hand, the prefactor A(t w ) is 
basically independent of the fitting-window. Another tricky point is the error-estimate 
for A(t w ) and a(t w ). It is clear that standard techniques (x 2 minimization) work poorly 
for tremendously correlated stochastic variables such as C(t, t w ) for successive t. In order 
to have a (hopefully) reasonable estimate we have turned to a Jack-Knife procedure on 
the fitted coefficients themselves. 

Our results for the prefactor A(t w ) and the exponent a(t w ) of the heat-bath 
dynamics in the D = 3 case, can be found in Figs. |U (T = 0.8), El (T — 0.6) and 
E1(T = 0.5). The results for t w < 8192 were obtained on a PC (see section^, while for 
t w > 8192 SUE data were used. Actually, we have found that the functional form © is 
not suitable for very small t w , thus there is a lower, temperature-dependent, cut-off on 
the t w values shown in Figs. IH El and El 

As far as the exponent a(t w ) is concerned (see the upper part of Figs. |U ElandE}, 
we have a significant growth for t w < 10 6 . For larger t w , a(t w ) is constant within 
errors. This (asymptotic?) value would correspond to the inverse of the large waiting 
time limit of Rieger et al. [9J aging exponent \(T,t w ). Interestingly enough, the values 
of a{t w ) for the largest achieved t w seems to be proportional to T c /T. On the other 
hand, the prefactor A(t w ) (see the lower part of Figs. |U El and EJ) clearly decreases in 
all the simulated t w range. A power law seems to be appropriate for this decay, the 
exponent being precisely the exponent x(T) in Eq.fjHl). Thus, we find x(0.7T c ) = 0.016, 
x(0.53T c ) = 0.0064, x(0.44T c ) = 0.0049. 

We have found similar results (see Fig. El) for the Metropolis dynamics of the 
Edwards-Anderson model with binary-couplings distributions in three dimensions (T = 
0.7T C ), four-dimensions (T = 0.6T C ). For all the simulated t w (< 10 6 ) we have found a 
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Figure 6. As in FigHfor 0.44T C 
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10' ! 



o O O 



10 c 



Figure 7. Top: exponent a(t w ) defined in Eq.|0 versus t w for the Metropolis 
dynamics of the Edwards- Anderson model in D = 3 at T = 0.8 (full circles) and 
D = 4 at T = 1.2 (empty circles). Bottom: same as in top panel for the prefactor 
A(t w ), defined in Eq.©. 



growing trend in the exponent a(t w ). As usual, the dynamics in four dimensions is faster 
than in three dimensions (notice how the a(t w ) is smaller but more rapidly growing). 
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Figure 8. Top: The spin correlation function as a function of t for several t w , obtained 
with the Metropolis dynamics of the 3D Edwards- Anderson model at T = 0.7T C (the 
lines are mere guide-to-eye). Data were obtained in a L = 33 lattice simulated on a PC. 
The plot concentrates on the region of large C in order to discuss the time-traslation 
invariant regime (see text). 



The prefactor A(t w ) decreases. A power-law seems to be adequate for this decay in 
D = 3, but this is less clear for the four- dimensional lattice. 

Up to now, we have only confirmed that the dynamics for t > t w , and t w > 10 6 
can be rather well described by the cross-over (201 IE] formula (j3J), and we have given 
an explicit form for the cross-over function <&(?/) (see Eq.©). We are thus predicting 
that C(t w ,t w ) should vanish in the infinite waiting-time limit: dynamic ultrametricity 
reduces to a trivial = statement in the t > t w time-sector. Therefore, non trivial 
statements about dynamical ultrametricity should regard times 1 <C t <C t w . In this 
time sector, the cross-over formula (j3J, predicts a time-translational invariant (fi = 0) 
power law decay, C(t,t w ) = $(0)/t 2 ^. This prediction is of course non-ultrametric. 

However, one could wonder about the presence of more than one time sector. In 
fact, as one can see in Fig. |H] for the Metropolis dynamics of the 3D Edwards- Anderson 
model at T = 0.7T C , small but measureable deviations from time-traslation invariance 
appear at t/t w ~ 10~ 3 . In order to explore the regime 1 t t w , we have introduced 
a subtracted correlation function: 

/ f \ !/«(*«>) 

C BuhB (t,t w ) = C(t,t w ) - A(t w ) 1 + — • (10) 
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Figure 9. Top: The subtracted correlation function defined in Ea. (jlOf) as a function 
of t/t w , for the Metropolis dynamics of the 3D Edwards- Anderson model at T = 0.7T C 
(the lines are mere guide-to-eye). The inset shows C su h s {t,t w ) as a function of time, 
for the different t w . Bottom: C su b s (t,t w ), from the same simulation as above, versus 
the (not dimensionless!) ratio t/t^ 6 . 



where A(t w ) and a(t w ) are of course the coefficients obtained in the fit to the functional 
form Q for each t w . To motivate it, let us recall that the quasi equilibrium regime 
is realized in the case where t w goes to infinity at fixed t. Naive scaling predicts that 
C(t + t w , t w ) goes to a non trivial function of t/t w when t w goes to infinity. However 
we have already remarked that we could have also a subaging contribution and find a 
non trivial behaviour in the region where t is of order with \i < 1. In this region we 
have 1 <C t <ti t w , but we may be not in the quasi-equilibrium regime as far as t goes 
to infinity and it is not fixed. One reason for studying C su b s (t, t w ) is that in the case 
of multitime sectors, the correlation function is the sum of contributions coming from 
each time sector and the subtraction help to identify the given time sector. 

We show in Fig.|H]Csubs(t, t w ) for D = 3, T = 0.8, as a function of t/t w (upper-part) 
and as a function of t/t^ 6 (lower panel). It is clear that C su b s can be described as a 
function of t/t w for t > 0At w , but that strong deviations are present for smaller times. 
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Figure 10. Double logarithmic plot of C(t,t w ) versus 1 + t/t w for the Metropolis 
dynamics of the infinite-dimensional z — 6 Viana-Bray model in (N = 10 6 , T — 0.8 = 
0.34T C ). The clearly measurable curvature indicates deviations from Eq.JHJ- 



The inset in the top panel of Fig. 01 shows as well that C sn b s (t,t w ) is not a function 
of t, as the cross-over formula (j3J) would suggest. On the other hand, the lower part 
of Fig. El shows that C S ubs(Muj) seems really a function" 1 " of t/t^ 6 during two decades 
(corresponding to the decay from C su bs = 0.18 to C su b s = 0.1). Up to our knowledge, this 
is the first time that a time sector different from /i = (time-translational invariant) 
and \i — 1 (full-aging), has been studied. A cautionary remark is in order, though, 
for one cannot exclude that the t/t^ scaling could stop to apply at much larger t w . 
Unfortunately, because of the small number of simulated samples, we have not been 
able to repeat this analysis with SUE data: the t/t w scaling and the t/t^ 6 scaling for 
Csubst look indistinguishable within error bars. 

It is clear that more work needs to be done in order to design an efficient protocol 
to study the different time-sectors. Yet, we hope that the reader will be convinced that 
it should be feasible. 

The Metropolis dynamics of the infinite-dimensional z = 6 Viana-Brey model is 
rather different. As we show in fig. the long-time decay of C(t,t w ) is not a power- 
law. Thus, it cannot be fitted with Eq. ([§]) . We have found (see fig. fTTj) that a fit 

+ Notice that the numerical value of the quotient t/t^ e is not invariant under a change of time units, so 
that one should really speak about t/(r^ A t^ 6 ), tq being the used time unit. Thus, we do not attribute 
any special significance to t/t^ 6 = 1. 
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Figure 11. Logarithm of of C(t,t w ) versus ^/log(l + t/t w )/ log t w (see Ea. (fTT|) ') 
for the Metropolis dynamics of the infinite-dimensional z = 6 Viana-Bray model in 
(N = 10 6 , T = 0.8 = 0.34T C ). The dashed horizontal line is the logarithm of the 
<?ea = 0.78 |2S]. 



to 



C(t,t w ) = S(t w ) exp 



-V(tv 



log (I + t/t v 
logt w 



is rather adequate. Moreover, and 77 (t^) seem to have a well defined limit for large 

t w , S(t w — > oo) being quite close to the Edwards-Anderson order-parameter q^A = 0.78 
that in this model has been computed [2E] at the one-step level of replica-symmetry 
breaking. If one believes that our results are almost asymptotic, so that the large t w 
correlation function truly is 



C(t, t h 



q EA exp 



-v 



log (! + */*„ 



log t v 



t > t, 



(12) 



the rather amazing conclusion is reached that this ultrametric system (from the point 
of view of statics [3]) is not ultrametric from the point of view of dynamics! It has been 
pointed out [TT] that dynamical ultrametricity implies (under reasonable hypothesis) 
statical ultrametricity, however our results suggest that dynamical ultrametricity is 
maybe not a necessary condition for the validity of RSB. The reader can check that the 
weaker property of separation of time scales holds: let t(C;t w ) be the time necessary 




Figure 12. Cu n k(t w , t w ) and qn n k(tw) versus t~ ' 08 , for the heat-bath dynamics of the 
3D Edwards- Anderson model at temperatures T = 0.7T C , 0.53T C and 0.44T C . Lines are 
linear fits. 



for the correlation function at waiting-time t w to reach the value C. One has 

i lim oo7pT4 = if Cl> ^- ( 13 ) 

It is somehow disappointing, though, that this property is expected to hold as well for 
the Langevin dynamics of the disordered ferromagnet jHj. 



4. The link overlap 



In recent years, a new picture of the low temperature phase of spin-glasses has been 
put forward [T3], the so-called TNT picture. These authors propose that the link- 
overlap defined in Eq. (|H|) should have a trivial distribution function (a Dirac delta) in 
the thermodynamic limit (for a recent study of the link-overlap and related quantities 
see Ref. |2Z|)- On the other hand, the spin-overlap defined in Eq.fjZJ), would have 
a non trivial distribution function, as predicted by RSB In contradiction with 
RSB, the TNT picture requires that the large waiting time limit of Cn^tw, t w ) be 
equal to the limit of qn-akitw) (since there would be only a possible value for this 
quantity!). This can be checked in our simulations. In Fig. [12] we show C\i n ^(t W) t w ) 
and qw-nkitw) obtained from SUE, for temperatures T = 0.7T C , T = 0.53T C and 
T = 0.44T C , as a function of t~ om . The extrapolation as a function of t~ ' 08 was 
suggested by the fact that giink(^) is roughly a linear function of this variable. Some 
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Figure 13. Top: Coefficients a(t w ) and b(t w ) (defined in Eq.lfHjl. versus t~ ' 17 , 
for the Metropolis dynamics of the 4D Edwards-Anderson model at temperatures 
T = 0.6T C . The fits have been performed in the window 0.1 < C 2 (t,t w ) < 0.3. 
Bottom: Cunk(t w ,t w ) and qunk(tw) versus fc^ ' 1 , for the Metropolis dynamics of the 
4D Edwards- Anderson model at temperatures T — 0.6T C . 



arguments for its validity has also been given in previous work [23 123 121] ■ We notice 
that q\ink(t w ) has a very mild temperature dependence. Furthermore, for large t w , 
Q]iak(t w ,T = 0.6) and qn n k(t w ,T = 0.5) are on the top of each other (within error 
bars), while q\mk(t w ,T = 0.8) is fast approaching them. Indeed, one would expect 
<?iink(^ = oo,T = 0.5) > qnnk{t w = oo,T = 0.6) > qn n k(tw = oo,T = 0.8), so that the 
three lines should cross. The fact that the infinite waiting-time limit of qn n k{t w , T = 0.5) 
is smaller than qn n k(t w = oo,T = 0.6) > q\i n k(t w = oo,T = 0.8), indicates that for 
T = 0.5 the data should probably extrapolate like t~ v with 77 < 0.08. On the other 
hand, C y m k(t m , t w ) shows a much stronger temperature dependence and, for T = 0.5 and 
T = 0.6, is basically t w independent. It seems plausible that, for T = 0.8, the limiting 
value of Ciink (t w , t w ) and qimk(t w ) will be fairly close. However, for T = 0.5 and 0.6 
the limits will be noticeably different unless the dynamics changes drastically at larger 
t w . Interestingly enough, in four dimensions and T = 0.6T C (see Fig. [T^J lower part), 
Cn n k(t w ,t w ) and q\mk{t w ) seems to be linear in t~ ' 17 , and to extrapolate to different 
values. 

Another interesting question regards the aging properties of the link correlation- 
function. From the droplet [Hj and TNT pictures of spin-glasses, one would not expect 
that Ciink (t, t w ) would age (at least, for large enough t w ). On the other hand, the RSB 
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Figure 14. C un k(*> t w ) versus C 2 (t,t w ), obtained in SUE for the heat-bath dynamics 
of the 3D Edwards- Anderson model at T = 0.7T C . The plot is restricted to the window 
of C 2 where a linear dependence is found for all the studied t w . The lines are linear 
fits. 

picture expects aging properties akin to the ones of C(t, t w ). To check for this, one could 
just look to Ciink (t, tw) as a function of C 2 (t,t w ) (see Fig. EJ). If one concentrates in a 
small window of C 2 , 0.2 < C 2 < 0.35 a linear description is perfectly adequate: 

CiinkO, t w ) = a(t w ) + b(t w )C 2 (t, t w ) . (14) 

The question now translates to the behavior of a(t w ) and b(t w ). The TNT and droplet 
pictures predict that b(t w ) tend to zero and a(t w ) tend to the large waiting time limit 
of Cii n k (t w , tw)- On the other hand, RSB predicts a non-vanishing limit of b(t w ). 

In Fig|T31 we show coefficients a(t w ) (bottom) and b{t w ) (top), as a function of 
t~ om , for the heat-bath dynamics of the 3D Edwards-Anderson model, at T = 0.5, 0.6 
and 0.8. Obviously, the lower T, the smaller is the range of t w shown, because C 2 (t,t w ) 
does not reach 0.2 for all the t w in our simulation time. The errors in a(t w ) and b(t w ) 
have been calculated with a Jack-Knife procedure. The coefficients seems to have a nice 
linear behavior in t~ ' 08 . One thus conclude that, unless a drastic change arise in the 
behavior of a(t w ) and b(t w ), the slope b(t w ) will not vanish asymptotically. This implies 
that Cjiok(tjt w ) should age as C 2 (t,t w ) does. One must acknowledge, however, that the 
conclusion is less sound for T = 0.8 than for the lower temperatures. In fact, for T = 0.8 
the infinite-volume extrapolation of a(t w ) (see bottom part of Fig. ITHjl is quite close to 
the large time-limit of Ci in k(t t0 , t w ), as TNT would predict. For the lower temperatures 
this extrapolations are clearly different. The data in four-dimensions support as well 
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Figure 15. Top: Coefficient b(t w ), defined in Ea. (|14(l obtained in SUE for the heat- 
bath dynamics of the 3D Edwards- Anderson model versus t~ ' 08 , for temperatures 
T = 0.7T C , 0.53T C and 0.44T C . Bottom: As in top part, for coefficient a(t w ), defined 
in Ea. (jl4|l . The lines are linear fits. We also plot in the vertical axis the infinite-limit 
extrapolation for C\i n ^(t w , t w ) (see Fig. ll2|l for all three temperatures. Notice that for 
T = 0.8 the infinite t w limit of a{t w ) and of Cn n y{t w ,t w ) seem compatible. 



the RSB prediction (see the top part of Fig[T3j) . 
5. Conclusions 

In this work, we report the results of a large scale Monte Carlo simulation of the heath- 
bath dynamics of the three dimensional Edwards- Anderson model with binary coupling 
at temperatures T = 0.7T C , 0.53T C and 0.44T C (see table [TJ. The long times achieved 
and the large lattices studied (60 3 ), have been made possible by the SUE machine. 
In addition, shorter but more precise* simulations are presented for the Metropolis 
dynamic of the same model in 3D, 4D and the infinite-dimensional Viana-Bray model. 

For the spin correlation-function in D = 3 and D = 4, we find that an slightly 
modified power law (see Eq.®), well describes the data for t > t w . This formula is 
identical to the cross-over like parametrization Q proposed by Rieger et al. EH- 
However, in the regime 1 <C t <C t w a different behavior is observed. A numerical 
procedure is proposed for the study of the different time-sectors believed to exist in 
spin-glasses dynamics |8j. Indeed, a time sector with characteristic exponent fi = 0.6 

* Due to the larger number of simulated samples. 
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is observed for the first time, we believe, although to firmly establish this result will 
require very precise simulations at still larger t w . In spite of this, the very slow t w 
evolution of A(t w ) and a(t w ) (see Eq.(JHJ)), makes us to believe that in thermoremanent 
magnetization measurements (usually restricted to the range lO -2 ^ < t < 10 2 t w ), a 
perfect full-aging behavior occurs, in agreement with a recent experiment ^Hj, and mild 
disagreement with older ones [H] (see also [IZj). Although it is somehow disappointing 
that the study of the (in)existence of more than two-time sectors in spin-glass dynamics 
should be restricted to simulations, we still think that it should be feasible. Nevertheless, 
one can always hope that new experimental techniques and protocols will be eventually 
able to explore this time-regime. 

In infinite-dimensions, or at least for the z = 6 Viana-Bray model, we have 
found a different scaling. The decay of the spin-correlation function is not a power 
law. Moreover, the limiting functional-form Eq. (jl2j) is not dynamically ultrametric. 
Although clearly more work is needed to establish this result, it suggests that dynamical 
ultrametricity is maybe not such an interesting property as previously thought. 

We have considered as well the aging properties of the link-overlap and the link 
correlation-function, both in three dimensions and in 4D. We have concluded that, unless 
a drastic change in the dynamical properties arises for t w larger than the here studied, 
the link correlation-function should age precisely in the same way as the spin correlation- 
function. This is in plain disagreement with the droplet and the TNT pictures of the 
spin-glass phase. However, one must acknowledge that the data at the highest studied 
temperature in three dimensions (T = 0.7T C ) are not incompatible with the TNT picture. 
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